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ABSTRACT 

Detailed observational characterization of transiting exoplanet systems has revealed that the spin- 
axes of massive (M > 1.2 M©) stars often exhibit substantial misalignments with respect to the orbits 
of the planets they host. Conversely, lower-mass stars tend to only have limited obliquities. A similar 
trend has recently emerged within the observational dataset of young stars’ magnetic field strengths: 
massive T-Tauri stars tend to have dipole fields that are ^10 times weaker than their less-massive 
counterparts. Here we show that the associated dependence of magnetic star-disk torques upon stel¬ 
lar mass naturally explains the observed spin-orbit misalignment trend, provided that misalignments 
are obtained within the disk-hosting phase. Magnetic torques act to realign the stellar spin-axes of 
lower-mass stars with the disk plane on a timescale significantly shorter than the typical disk lifetime, 
whereas the same effect operates on a much longer timescale for massive stars. Cumulatively, our 
results point to a primordial excitation of extrasolar spin-orbit misalignment, signalling consistency 
with disk-driven migration as the dominant transport mechanism for short-period planets. Further¬ 
more, we predict that spin-orbit misalignments in systems where close-in planets show signatures of 
dynamical, post-nebular emplacement will not follow the observed correlation with stellar mass. 


1. INTRODUCTION 

Of all ideas in planetary science, few have stood the 
test of time better than the “Nebular Hypothesis”, origi¬ 
nally proposed in the 18th century by Kant and Laplace 
(Kant 1755; Laplace 1796). The impetus for such a model 
was the near-aligned configuration of the planetary orbits 
with each other and with the Sun’s spin axis; the net 
angular momentum of the the planets differs in direction 
from the Sun’s spin by about 7 deg^ (Beck & Giles 2005). 
Centuries of refinement (with help from astronomical ob¬ 
servations) have resulted in the modern picture whereby 
a dense molecular cloud core collapses under its own self 
gravity (Shu et al. 1987; McKee & Ostriker 2007) to form 
a star encircled by a disk of gas and dust. 

The earliest descriptions of molecular core collapse 
were naturally the most simplistic, supposing that the 
star and its disk both inherit similar angular momen¬ 
tum directions. More recent work (Goodman et al. 1993; 
Caselli et al. 2002; Bate et al. 2010) has added a layer of 
complexity to the story by noting that turbulence in col¬ 
lapsing cores implies that the last bits of material falling 
onto the disk do not necessarily share the same angular 
momentum direction as the star. Despite this apparent 
tendency towards slight misalignment, the mutual grav¬ 
itational torque between star and disk is likely strong 
enough during the earliest stages to stave off any signifi¬ 
cant star-disk misalignment arising from core turbulence 
(Spalding et al. 2014). 

Long before the first exoplanetary detections (Mayor 
& Queloz 1995; Marcy & Butler 1996), the theory behind 
protoplanetary disks was already fairly mature (e.g., Gol- 
dreich & Tremaine 1979, 1980; Lin & Papaloizou 1986; 
Ward 1986). An early prediction from this field was that 
angular momentum exchange between disks and their 
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^ Of course, 7 deg is very different from zero and so is still in 
need of an explanation. 


embedded planets should give rise to planetary migra¬ 
tion towards shorter-period orbits. More recent work 
has subjected this idea to extensive numerical (Nelson 
et al. 2000; Rice et al. 2008) and analytic (Tanaka et 
al. 2002) analyses, with the general picture of migration 
holding up as an expected outcome (for a recent review, 
see Kley & Nelson 2012). Sure enough, the earliest days 
of exoplanet-hunting revealed a considerable population 
of hot Jupiters, planets with about the mass of Jupiter 
but with orbital periods of a few days. Conventional 
planet formation theory (Pollack et al. 1996) suggests 
that these planets must form at several AU, beyond the 
snow line, and must have migrated towards their present- 
day orbits. Accordingly, disk-driven migration seemed 
an attractive mechanism by which giant planets may be 
delivered to close-in orbits. 

Until recently, it was impossible to tell whether these 
systems, supposed to migrate through a planar disk, were 
aligned or misaligned with their stars. This property, re¬ 
ferred to variably as obliquity or spin-orbit misalignment, 
has now fallen within the observational capabilities of ex¬ 
oplanetary astronomy by way of the Rossiter-McLaughin 
effect (Rossiter 1924; McLaughlin 1924; Winn et al. 
2005). Currently, measurements of misalignments are 
most common for the orbits of hot Jupiters, wherein the 
findings have revealed that a substantial fraction of such 
planets follow orbits with significant obliquities (Winn 
et al. 2010; Albrecht et al. 2012). Indeed, misalign¬ 
ments range all the way from prograde aligned to ret¬ 
rograde anti-aligned. However, the degree of misalign¬ 
ment exhibits a clear dependence on stellar mass, with 
the most extreme, retrograde (circular) orbits only ap¬ 
pearing around stars with mass greater than M > 1.2 Mq 
(Figure 1). 

In light of the existence of significant obliquities, many 
authors (e.g., Albrecht et al. 2012; Lai 2012; Dawson 
2014; Storch et al. 2014; Petrovich 2015) have sought al¬ 
ternatives to disk-driven migration as a production mech- 
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Fig. 1. — The observed projected angle between the stellar spin 
axis and orbital plane of circular planetary orbits (e < 0.1; red, 
solid points) and eccentric orbits (e > 0.1; black, faint points) for 
stars of given masses. There exists a clear distinction between 
low-mass stars (M < 1.2 Mq), which display moderate-to small 
misalignments (especially among circular systems), and the more 
massive stars (fif > 1.2 Mq), which exhibit misalignments ranging 
all the way from retrograde-aligned to prograde-aligned. Measure¬ 
ments of magnetic field strengths (Gregory et al. 2012) have re¬ 
vealed that, among T-Tauri stars, lower-mass stars (similarly, cor¬ 
responding to M < 1.2 — 1.4 M0) possess a much stronger surface 
dipole field than do their higher-mass counterparts. Specifically, 
low-mass stars possess fields of ~1 kGauss in contrast to more mod¬ 
est ~0.1kGauss for higher-mass T-Tauri stars. The misalignments 
data were obtained from exoplanets. org and follow the discussion 
of Albrecht et al. (2012). 

anism for inclined hot Jupiters. Specifically, in contrast 
to the ‘smooth’ picture of disk-driven migration, there 
exists a separate class of migration mechanisms which 
occur after the dissipation of the protoplanetary disk. 
Post-disk migration must invoke dynamical interactions 
to excite near-unity eccentricities, bringing the perias- 
tron close enough to the central star such that tidal 
forces may take over and circularize the orbit (Ford & Ra- 
sio 2006). The necessary perturbations are hypothesized 
to arise from processes such as planet-planet scattering 
(Ford & Rasio 2008; Nagasawa et al. 2008; Beauge & 
Nesvorny 2012), Kozai resonance with a perturbing com¬ 
panion (Wu & Murray 2001; Fabrycky & Tremaine 2007; 
Naoz et al. 2011) or secular chaotic excursions (Lithwick 
& Wu 2012). 

Despite the natural tendency of dynamical interactions 
to excite inclinations, their occurrence rate appears to 
be insufficient to explain the current data (Dawson et al. 
2015). Conversely, disk-driven migration is expected to 
be almost ubiquitous but has long been assumed to give 
rise to low-obliquity systems because of the tendency for 
an isolated collapsing core to form an aligned star-disk 
system. Flowever, stars do not form in isolation. On 
the contrary, binary fraction probably lies somewhere be¬ 
tween half and unity during the early, disk-hosting stage 
of stellar evolution (Chez et al. 1993; Kraus et al. 2011; 
Marks & Kroupa 2012; Duchene & Kraus 2013). The 
potential for neighbouring stars to influence disk orienta¬ 
tion has actually been recognized for some time (e.g. Lar¬ 
wood et al. 1996), but it was first suggested as an expla¬ 
nation for spin-orbit misalignments by Batygin (2012). 
In this scenario, the gravitational torque exerted by a 
distant companion star causes the disk to process about 


the plane of the binary system. Observations suggest 
that binary planes do not correlate with disk orientation 
(Stapelfeldt et al. 1998; Koresko 1998; Jensen & Akeson 
2014) and so the disk-precession tends to torque the star 
and disk out of alignment. 

We recently advanced the so-called “disk-torquing” 
framework by taking into account the gravitational 
torques communicated between the star and its disk 
(Spalding & Batygin 2014). Specifically, the star tends 
to process about the disk’s plane at a rate that decreases 
with time owing both to stellar contraction and disk 
mass-loss. In the earliest stages, the star is well enough 
coupled to the disk to adiabatically trail the disk as it 
processes. However, once the precession rate of the star 
about the disk is roughly commensurate with that of the 
disk about the binary plane, a secular resonance is en¬ 
countered which is capable of exciting star-disk misalign¬ 
ments occupying the entire observed range. Such a pic¬ 
ture has been independently corroborated in other work 
and appears to be a robust result across various levels of 
approximation (Batygin & Adams 2013; Lai 2014; Xiang- 
Gruess & Papaloizou 2014). 

From the above discussion, it appears that spin-orbit 
misalignments may arise naturally out of either class of 
proposed migrationary scenario. What is less obvious is 
how to reconcile the migration pathway with the mass- 
dependence of the obliquity distribution (Figure 1). Pre¬ 
vious attempts to explain the trend (Winn et al. 2010; 
Lai 2012) have largely focussed on the violent migration 
pathway. Specifically, it has been proposed that dynam¬ 
ical encounters misalign orbits, after which, strong tidal 
coupling between low-mass stars and their planets (ow¬ 
ing to an extended convective region) re-aligns the orbits. 
These hypotheses rest upon a number of assumptions, 
one of which being that the tidal damping rate of obliq¬ 
uity greatly exceeds that of orbital semi-major axis (Lai 
2012). Additionally, the presence of misaligned multi- 
transiting systems (Huber et al. 2013) has yet to be given 
a viable explanation within the violent framework. 

The disk-torquing mechanism is expected to misalign 
all planets in the system by a similar amount, at least in 
the inner disk. However, to date, the mass-misalignment 
trend has not been given an adequate explanation. Here, 
we provide such an explanation by utilizing an addi¬ 
tional piece of evidence in the form of the magnetic field 
strengths of young, disk-hosting stars. Specifically, re¬ 
cent observations (Gregory et al. 2012) have revealed 
that low-mass T-Tauri stars possess dipole field strengths 
of order ~ 1 kG, an order of magnitude greater than 
their higher-mass counterparts (^ 0.1 kG). This trend 
appears to continue into the more massive Herbig AeBe 
class of star-disk systems (M ~ 1.5 — 8 M 0 ) (Alecian 
et al. 2013), suggesting the dipole strength to be a ro¬ 
bust mass-dependent feature, intrinsically connected to 
PMS evolution. Grucially, the transition from weaker 
to stronger fields occurs at a similar mass to that defin¬ 
ing the mass-misalignment trend (~ 1.2M 0 ). In this 
work, we suggest that the two trends are causally linked: 
stronger magnetic fields of low-mass stars erase primor¬ 
dial misalignments that their higher-mass, weaker-field 
counterparts retain. 

In order to test this idea, we must hrst quantify the 
torques communicated between a tilted dipole and a disk. 
We are not the first to suggest that magnetic torques 
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might influence misalignments (Lai et al. 2011) but our 
goal is to couple the star-disk magnetic torques to the full 
disk-torquing picture, i.e., a freely precessing star-disk 
system. Consequently, we develop a more complete pic¬ 
ture of the various torques involved using semi-analytic 
methods, whilst taking advantage of the conclusions of 
various numerical studies (mostly on aligned disks, e.g. 
Ghosh & Lamb 1978; Armitage & Clarke 1996; Uzdensky 
et al. 2002). Highly sophisticated numerical models exist 
for the study of star-disk interactions (e.g. Romanova 
et al. 2012), however, detailed simulations must be in¬ 
tegrated in full throughout each stellar rotation period. 
Even at the current state-of-the-art, full 3D calculations 
can only be carried out within a reasonable length of 
computer time for ^ 1000 orbital timescales (depend¬ 
ing upon the precision). In contrast, we seek to model 
the global evolution of the star-disk system, including 
time-dependent disk mass and stellar contraction, over 
multi-Myr timescales. Considering also that current ob¬ 
servations only constrain the topologies of T-Tauri field 
strengths to within an order of magnitude or so, we shall 
construct our model such that the level of detail is com¬ 
mensurate with that of the observations. 

In this paper, we begin with a description for the time 
evolution of the star and disk. Next we provide a re¬ 
hashing of the purely gravitational torques as derived 
in Spalding & Batygin (2014), derive the various mag¬ 
netic torques communicated between the disk and star, 
and subsequently present the results obtained through 
numerical integration of the star-disk system, complete 
with gravitational and magnetic torques. We conclude by 
discussing the implications of the result for the acquisi¬ 
tion of spin-orbit misalignments within the disk-torquing 
framework. 


2. MODEL 

Our goal is to describe the spin-axis dynamics of a 
star, possessing a dipolar held that is encircled by a pro¬ 
toplanetary disk. We suppose the star-disk system to be 
orbited by a companion star. In what follows, we derive 
the analytical forms of the various torques inherent to the 
system. Throughout the entire calculation, we assume a 
hierarchical configuration. Specifically, we assume that 
the central star and binary companion do not influence 
each other directly. Rather, the companion torques the 
disk (gravitationally), which in turn interacts with the 
central star (gravitationally and magnetically). In order 
to improve clarity, we adopt the following convention for 
variables with identical symbols: quantities referring to 
the disk are primed, those referring to the central star 
are marked with a tilde, and those referring to the com¬ 
panion star are given an over-bar (e.g. semi-major axes 
will be defined as a', 5, and d respectively). 

The binary companion is prescribed a mass M and 
an orbit with semi-major axis d which is much greater 
than the outer disk edge « 30 AU. The resulting 
gravitational torques acting upon the disk cause it to 
process at a frequency that depends upon the mass of 
the companion, its semi-major axis and its inclination 
relative to the disk-plane (Spalding & Batygin 2014). In 
the absence of star-disk interactions, the precession of 
the disk would simply cause it to tilt out of alignment 
with its host star. However, the central star can inter¬ 
act with its disk by way of several physical processes. 


Specifically, young stars rotate rapidly enough to pos¬ 
sess a considerable centrifugal bulge at their equators. 
This bulge allows for gravitational coupling between the 
disk and star, the dynamics of which having been derived 
in detail elsewhere (Batygin & Adams 2013; Spalding & 
Batygin 2014). In addition, observations of young stars 
have revealed the presence of significant magnetic fields 
(Johns-Krull 2007; Gregory et al. 2012) which facilitate 
angular momentum transfer between star and disk in ad¬ 
dition to the gravitational influences. 

Our ultimate goal is to develop a theoretical framework 
by which we may test the hypothesis that differences in 
stellar magnetic field between high and low-mass stars is 
the dominant driver of the observed mass-misalignment 
trend in the current Rossiter-McLaughlin dataset (Winn 
et al. 2010; Albrecht et al. 2012). We follow a semi- 
analytic, parameterized framework in deriving the rele¬ 
vant equations. 

As the system ages, the star contracts and the disk 
loses mass. Accordingly, before calculating the star-disk 
torques, we must first provide formulations of the phys¬ 
ical evolution of the disk and the central star. Within 
such a framework we then describe our calculations of 
magnetically-facilitated tilting of the stellar-spin axis and 
star-disk gravitational coupling. Torques arising from ac¬ 
cretion are neglected in this work because their influence 
within a similar physical framework has been examined 
elsewhere (Batygin & Adams 2013) and found to be in¬ 
significant. 


2 .1. Physical Evolution of the Protoplanetary Disk and 
the Stellar Interior 

Typically quoted lifetimes of protoplanetary disks fall 
in the range ^1 — 10 Myr, with some recent evidence 
in support of longer-lived systems (Haisch et al. 2001; 
Williams & Cieza 2011; Bell et al. 2013). We parameter¬ 
ize disk mass evolution as (Laughlin et al. 2004): 


Aldisk 


<sk 

1 4“ t/rdisk 


( 1 ) 


The time derivative of Mdisk approximately represents 
the accretionary flow. Following Spalding & Batygin 
(2014), we note that observations (Hartmann et al. 2008; 
Herczeg & Hillenbrand 2008; Hillenbrand 2008) are best 
matched by adopting an initial disk mass, = 5 x 

10“^ Mq and dissipation timescale Tdisk = 5 x 10“^ Myr. 

For simplicity, we approximate the central star with 
a polytrope of index ^ = 3/2 (appropriate for a fully 
convective object; Chandrasekar 1939). A polytropic 
body of this index possesses a specific moment of iner¬ 
tia I = 0.21 and a Love number (twice the apsidal mo¬ 
tion constant) of ^2 = 0.14. In contracting along their 
respective Hyashi Tracks, T-Tauri stars derive most of 
their luminosity from the release of gravitational poten¬ 
tial energy. We describe the radiative loss of such binding 
energy as (Hansen & Kawaler 1994): 




GMf dR^ 
2Rl ~dl' 


(2) 
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with a solution, 


R* = (R’l) 



24TraT^s 

GM^{R0)3 



( 3 ) 


For definiteness, we match the numerical evolutionary 
track of an M* = 1 Mq star (Siess et al. 2000) by as¬ 
suming an initial radius of i?2 — 4i?Q and an effective 
temperature of Tgs = 4100 K. 


2.2. Disk-Binary Gravitational Interactions 

The binary companion interacts with the disk through 
gravitational torques alone. To capture the long-term 
behaviour of the angular momentum exchange, we uti¬ 
lize the secular approximation (Murray & Dermott 2000; 
Morbidelli 2002) whereby the torques exerted by the 
companion on the disk are equivalent to those communi¬ 
cated by a massive wire sharing the companion’s orbital 
elements. As a consequence of self-gravity and hydrody¬ 
namic pressure forces, the disk retains coherence under 
the influence of such torques, acting as a rigid body (Lar¬ 
wood et al. 1996; Batygin et al. 2011; Xiang-Gruess & 
Papaloizou 2014). Additionally, we neglect any dynami¬ 
cal influence of disk-warping. 

We adopt a Hamiltonian framework in describing the 
dynamics. The gravitational torques acting between a 
star, disk and companion have already been comprehen¬ 
sively derived and discussed elsewhere (Spalding & Baty¬ 
gin 2014). Accordingly, here we provide only a brief out¬ 
line of the corresponding computation. 

To begin, we introduce the scaled Poincare action- 
angle coordinates, defined in terms of the disk-star mu¬ 
tual inclination f}' and the disk’s argument of ascending 
node Vi', 


Z' = 1- cos(/30 z' = -n'. (4) 


the pre-main sequence. Of course, these complications 
will lead not only to a time-dependent v, but also to sec¬ 
ular variations in Z'. This point is not crucial to the 
problem at hand, so we leave analysis of such dynamic 
processes to future work and maintain a constant v and 
Z' throughout. 


2.3. Disk-Star Interactions: Gravity 

Having prescribed the secular evolution of the disk ow¬ 
ing to the companion (constant, rigid-body precession), 
we now describe the processes by which torques are com¬ 
municated between the star and the disk. Observations 
of T-Tauri stars have revealed that they spin with periods 
ranging between ~ 3 — 10 days (Herbst et al. 2007; Lit- 
tlefair et al. 2010; Bouvier 2013). Such high spin rates 
(^within a factor of ten of break-up rotation) lead to 
the development of a significant centrifugal bulge on the 
stellar equator. When the star and disk are misaligned, 
this bulge results in gravitational torques that force a 
precession of the stellar spin pole about the disk plane 
(analogous to a top spinning on a planar table). 

The dynamics of a spheroidal star and the gravitational 
influence arising from its rotationally-derived equatorial 
bulge may be approximated to high precision by con¬ 
sidering the inertially-equivalent picture of a wire with 
mass m in circular orbit with semi-major axis a around 
a point mass. Respectively, the appropriate mass and 
semi-major axis are given by (Batygin & Adams 2013) 


m = 


AGk2 



a = 


_ 9/2 GM* 



( 7 ) 


The appropriate Hamiltonian describing the compan¬ 
ion’s influence upon the disk is then 


3«/ut 

M 

( <ut ^ ^ 

r 

8 


1 d ) 

2 


where riout is the mean motion of the gas at the outer 
edge of the disk. Crucially, U contains no explicit de¬ 
pendence on z'. Therefore, the disk-binary inclination 
is a constant of motion. Appropriately, in our analysis, 
we carry out our simulations within a reference frame 
co-precessing with the disk, at a rate v = dz'/dt, given 
by 


V = 


dU 

IdZ' 


8 M* 


a'out 

a 


3 


1-Z' 


( 6 ) 


Boosting into such a frame is equivalent to subtracting 
vt from the argument of ascending node of the disk. 

Owing to the arbitrary nature of choosing the com¬ 
panion’s orbital parameters, we prescribe v and Z' in 
our problem independently. It is worth noting that the 
picture whereby a single companion remains on a cir¬ 
cular orbit throughout the entire disk lifetime is highly 
idealized. In reality there may exist multiple companions 
and/or companions might be gained and lost throughout 


With this prescription, the standard perturbation tech¬ 
niques of celestial mechanics can be applied to the 
spheroidal star (Murray & Dermott 2000; Morbidelli 
2002 ). 

By working in a frame co-precessing with the disk, 
we introduce a time-dependence to the Hamiltonian re¬ 
sulting in an apparent linear increase in the disk’s argu¬ 
ment of ascending node of v t. The Hamiltonian can be 
made autonomous by employing a canonical transforma¬ 
tion arising from the following generating function of the 
second kind (Goldstein 1950): 


<G ,2 = {z- vt) $, (8) 

where (j) = {z — vt) is the new angle and the new momen¬ 
tum is related to the old one through: 


~ _ dG2 
~ hi 


= 4). 


(9) 


In addition to removing explicit time dependence, we 
scale the Hamiltonian by v. Following the transforma¬ 
tions described above, the Hamiltonian takes on the fol¬ 
lowing form: 
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•H = -$+— 3($-2)$ 

-3(2 + 3($-2))$cos2(/3') 

+ 6 sin(2/3') ($ — l) a/( 2 — $)$ cos(0) 

+ 3 sin^(/?')($ —2)$ cos(2(/)) . 


R = ( — 


Po^(cos(0)) cos(/3') 


— sin(/3') I sin('i/’) sin(i^) 


( 10 ) 


■ cos('ip)cos((j)) )p^^(cos(0)) 


where, 


(13) 


where is the equatorial stellar surface field and P™ are 
associated Legendre polynomials. Within such a frame¬ 
work, the tilted dipole components are as follows, 


A= ^ ^ in -^disk Qjn \ /i 1 \ 

- 8 ^ ^ 

and is the mean motion at the disk’s inner edge. 

The purely gravitational dynamics described by the 
above Hamiltonian (10), together with the physical evo¬ 
lution of the star and disk, give rise to the excitement of 
mutual inclination between the central star and its disk 
(Spalding & Batygin 2014). However, additional physi¬ 
cal mechanisms must exist in order to explain the mass- 
misalignment trend. The main hypothesis of our paper 
is that the dominant driver of such a trend is the mass- 
dependence of T-Tauri dipole held strengths. Accord¬ 
ingly, next we present our derivations of the magnetic 
torques and the differential equations used in determin¬ 
ing their inhuence on stellar spin dynamics. 


2.4. Magnetic Torques 

In order to prescribe the magnetic disk-star interac¬ 
tions, we consider a T-Tauri star possessing a purely 
dipole magnetic held, whose north pole is aligned with 
the stellar spin axis. A pure dipole is modeled because 
the octupole component in real systems falls oh much 
faster with distance (oc 1/r®) than the dipole (oc 1/r^). 
Accordingly, at the position of the inner edge of the disk 
(^ 10 i?*), the octupole components have been attenu¬ 
ated to a factor of ^ 10“^ relative to the stellar surface 
held (r = ii*) as opposed to the ^ 10“^ attenuation suf¬ 
fered by the dipole component. Observations constrain 
the surface octupole held to diher from the dipole held 
by little more than a factor of 10 (higher or lower) and 
so the octupole component is almost always negligible 
when considering disk-star torques. The stronger dipoles 
of low-mass stars directly lead to a greater magnetic in¬ 
teraction with their disks than for high-mass stars. 

In the region of interest (i.e. in the domain of the disk), 
the stellar held is current-free and can be expressed as 
the gradient of a scalar potential: 



-I- cos{(j) — ij) sin(/3') sin(0)) 


(cos/3'sin(6») 

— cos((j) — V’) sin(/3') cos(0)) 

= {sin{l3’)sm((j) - Ip)). (14) 

which describe the vector held 

Bdip — B.j.Gj. BqGq -j- BfpGfj). (1^) 

If we assume that the disk material is in Keplerian or¬ 
bit about the star and for now suppose the star and disk 
are aligned, there exists a well-known expression for the 
corotation radius a(.o|/3'=o = (GM*/a;*)3, the radius at 
which relative angular velocity between the stellar mag¬ 
netosphere (w*) and the disk material is zero. Now con¬ 
sider tilting the star with respect to the disk. The az¬ 
imuthal motion of the stellar magnetosphere at the disk 
plane becomes reduced in such a way that the corotation 
radius is modulated as follows 


/ GM^ y 
Va;2 cos2(^')y 


(16) 


In other words, as the star tilts, the proportion of the disk 
which is effectively super-rotating relative to stellar rota¬ 
tion increases. At radii greater or smaller than the coro¬ 
tation radius, Keplerian shear results in relative motion 
between the stellar magnetosphere and the fluid compris¬ 
ing the disk (Figure 2). 

Owing to thermal ionization of alkali metals within the 
inner regions of the disk, motion of disk material relative 
to the stellar magnetosphere ‘drags’ the held lines, in¬ 
ducing magnetic helds. The evolution of induced helds 
is governed by the induction equation, 


Bdip = -VK 


( 12 ) 


dt 


(17) 


To retain generality, we take the held to be tilted at 
an angle pi' with respect to a spherical coordinate system 
(r , 6 , (p) into a direction specihed by an azimuthal angle 
fj: 


where B is the total magnetic held, i] is the magnetic 
diffusivity and 


V = Uif 60 — a; X r 


(18) 
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is the motion of the disk fluid in the frame rotating with 
the stellar magnetosphere^ i.e., we imagine the stellar 
magnetosphere to be held fixed, with the disk moving 
within it. The Keplerian velocity vk is azimuthally- 
directed whereas stellar rotation, arising from stellar an¬ 
gular velocity lo has both azimuthal and vertical compo¬ 
nents. 

A full, time-dependent solution of equation (17) is com¬ 
putationally difficult and so we adopt an approximate, 
semi-analytic approach. Specifically, we simplify the pic¬ 
ture by noting that the disk fluid velocity v possesses two 
separate components in the rotating frame. The first is 
an azimuthal component, arising from motion of disk ma¬ 
terial in Keplerian orbit relative to the azimuthal mag- 
netospheric motion from stellar rotation. Second, there 
is a vertical component, where ‘vertical’ refers to normal 
to the disk’s plane. Vertical motion occurs owing to the 
tilt of the star with respect to its disk, which causes the 
stellar magnetosphere to be dragged vertically through 
the disk every stellar rotation period (best imagined for 
a star tilted by 90deg; red diagram in Figure 2). In a 
frame rotating with the star, this vertical dragging ap¬ 
pears as a vertical motion of disk material (a vertical 
component to v) and therefore constitutes an additional 
source of magnetic induction (see section 2.4.2). 

For a small region on either side of the corotation ra¬ 
dius, the relative azimuthal motion may be sufficiently 
small to allow a steady state to exist between magnetic 
field dragging and slippage (Matt & Pudritz 2004) . How¬ 
ever, everywhere outside of this small region (of annular 
thickness ^ 0.01 a(.o), the diffusive timescale is longer 
than the dragging timescale, leading to the unbounded 
inflation of magnetic field lines. Such inflation can¬ 
not physically continue indefinitely and so some mech¬ 
anism must act to dissipate magnetic energy on orbital 
timescales. 

Analytic and numerical models of the aligned star-disk 
configuration find that the most likely dissipative pro¬ 
cess is magnetic reconnection, whereby the magnetic field 
lines are stretched azimuthally until the induced field is 
of the same order as the vertical stellar field. At this 
point the magnetic field lines break and reform (Livio 
& Pringle 1992; Uzdensky et al. 2002). Such a process 
is intrinsically non-steady, however. Owing to the short 
timescale of reconnection (^ the orbital period) com¬ 
pared to disk evolution, we may average over each orbit 
and consider a steady magnetic torque to result between 
the star and disk. 

We shall proceed by considering separately the fields 
generated from azimuthal relative motion and those orig¬ 
inating from vertical (in the disk’s frame) relative mag- 
netospheric motion. The form of equation (17) ensures 
that the fields induced by vertical and azimuthal veloci¬ 
ties back-react upon each other (through the v x B term) 
and so such separation of vertical and azimuthal induc¬ 
tion is not strictly accurate. However, as we are seeking 
an approximate model for the magnetic star-disk inter¬ 
actions, we proceed with the picture whereby the two 
components act independently. 


2.4.1. Azimuthal Induction 


Relative azimuthal disk motion induces an azimuthal 
field^ through flux-freezing (Armitage & Clarke 1996; 
Agapitou & Papaloizou 2000). This induced field, is 
represented as a fraction (also called a pitch angle), 7 of 
the component of the stellar dipole field perpendicular to 
the disk’s surface. In the case of a thin disk in spherical 
coordinates, the ‘vertical’ field is well approximated by 
the (negative) 0—component of the stellar field, Bg, at 
the disk mid plane (0 = 7 ’‘/ 2 ). Thus, = jBg, where 
the subscript ‘i’ refers to ‘induced’. 

As mentioned above, 7 is unable to instantaneously 
greatly exceed unity, but rather, magnetic reconnection 
reduces 7 to 7 ^ 0 each orbital period, allowing the 
field to be re-wound. We average over each reconnec¬ 
tion timescale and consider the star-disk torque to act 
equivalently to a steady torque of azimuthal pitch angle 
7 ~ 1. The force per unit volume is given by the Lorentz 
equation 


f = JxB. (19) 

In the case of azimuthal induction, the current density 
J = (l//ro)V X Bind may be considered to arise from the 
variation in magnetic field going from outside the disk 
vertically into the plane of the disk over some length 
scale S. In such a case, the current induced is radial and 
within the plane of the disk. Effectively, this statement 
is equivalent to setting the V operator equal to a vector 
of magnitude 1/5 in the vertical direction. Assuming 
5 ^ h where h is the disk scale height, the final torque 
per unit area on the disk is obtained approximately by 
multiplying equation (19) by <5, effectively integrating the 
torque over the vertical dimension. 

Angular momentum transport among neighbouring an¬ 
nuli of the disk is facilitated by the propagation of bend¬ 
ing waves (Foucart & Lai 2011), disk viscosity (Larwood 
et al. 1996) and through disk self-gravity (Batygin et 
al. 2011 ), all of which generally occur over a shorter 
timescale than does stellar tilting. As a result, the effec¬ 
tive moment of inertia of the disk in response to magnetic 
torques is much greater than that of the star, as was the 
case above in the gravitational picture. Accordingly, in 
the calculations which follow, we consider only the back- 
reaction of the torques, in determining the dynamics of 
the star (which effectively introduces a minus sign), with 
the disk’s dynamics being forced solely by the perturbing 
companion. 

From considering only the azimuthal field induction, 
arising from the penetration of stellar flux into the upper 
and lower surfaces of the disk, a torque per unit area of 


T = —V X 


'jBgej. X ^dip\g^^i2 


( 20 ) 


is communicated between the disk and the star. It is 
appropriate to display these torques in terms of their 
Cartesian components in the disk frame (i.e. points 
along the disk angular momentum axis) before integrat¬ 
ing over the entire disk. We find that the x, y and z 


^ Physically, the motion of the disk material relative to the 
background field generates a current within the disk. For purely- 
azimuthal fluid motion, the current is radial. This radial current in 
turn induces its own an azimuthal magnetic field (in the ideal MHD 
case), which leads the background field lines to appear stretched. 
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torques arising from equation (20) within the region in¬ 
terior to corotation (a < Oco) are given by, 


Tx 


Tz 


2 7rr (ago 
3 al 

2 7rr (ago 
3 al 

AttT 





sin(/3) cos(/3) cos('!/') 
sin(/3) cos(/3) sin(?/;) 
cos(/3)^ 


( 21 ) 


Once again, converting to Cartesian components and 
integrating, we obtain the torques arising from radial flux 
penetration in the form, 


Tx=0 

Ty=0 

SiraT . ,„,2 
Tz = - 3 — sm(/3) 


(25) 


where the variable T, a measure of stellar magnetic mo¬ 
ment, is defined for ease of writing as 

td2 jd6 

T=^^. ( 22 ) 

Mo 

Analogous torques generated by the region exterior to 
corotation are similar in functional form but instead of 
integrating the torques between Oin and Oco, we integrate 
outwards from Oco and note that the outer edge of the 
disk is sufficiently large as to be approximately equivalent 
to integration out to a' —)■ oo. Accordingly, the torques 
arising from outside of corotation are given by 


2 TT 'T' 1 

— -sin(/3) cos(/3) cos(V') 

^ ^co 

2 TT 1 

— -^ sin(/3) cos(/3) sin(V') 

^ ^CO 

(23) 

“co 

It is important to notice that the above torques are 
null for a star-disk inclination of /3 = 7r/2. The torque 
vanishes in such a scenario because in the picture con¬ 
sidered thus far, at Bg = 0 no flux penetrates the upper 
and lower surfaces of the disk. In other words, the null 
torque arises as an artifact of assuming a razor-thin disk. 
In reality, the disk possesses a finite scale height h. For 
small /3, most of the stellar magnetospheric flux pene¬ 
trates the upper and lower surfaces of the disk, validating 
the razor-thin model. However, as /3 is increased, more 
stellar flux penetrates the inner edge of the disk. This 
flux penetration leads to additional azimuthal twisting 
of field lines within a small annular region at the disk 
truncation radius (Figure 2). 

Whereas the induction considered previously generated 
a radial current, the induction owing to horizontal pen¬ 
etration of flux in the inner disk wall causes the field 
strength to vary along the radial direction. By the curl 
operator in equation (17), radial variation loosely trans¬ 
lates to the generation of a vertical current. If we again 
suppose that the induced held is of a similar magnitude 
to the background stellar held, the current generated 
is that given by a change in magnetic held strength of 
Br(r = Oin) (the radial stellar held evaluated at the trun¬ 
cation radius) over a radial length scale <5^- Once again, 
we integrate over this small region, assuming Sr/a-m 1, 
to hnd a torque per unit area (where the area is now the 
inner face of the disk) of 



1 

Mo 


Br&g X 


(24) 


where a = hmja-m is the aspect ratio evaluated at the in¬ 
ner edge of the disk. Its value is likely to be about ~ 0.1 
(Armitage 2011), although it may indeed be slightly 
larger at the inner disk edge as a result of thermal ex¬ 
pansion from strong stellar irradiation. For dehniteness, 
we set a = 0.1 in our calculations. 

2.4.2. Vertical Induction 

Many previous works have considered the torques aris¬ 
ing from azimuthal held dragging within accretion disks 
(Ghosh & Lamb 1978; Livio & Pringle 1992; Armitage & 
Clarke 1996; Agapitou & Papaloizou 2000; Uzdensky et 
al. 2002; Matt & Pudritz 2004; Lai et al. 2011). However, 
a feature which has been omitted from previous works is 
the fact that if a stellar magnetosphere is rotating at an 
inclination relative to its disk, there exists a component 
of relative star-disk motion which forces the held lines to 
be dragged vertically through the disk. This process is 
best imagined in the case of a star inclined by tt/2 to a 
disk, i.e., spinning on its side. In such a case, consider 
the situation in a frame co-rotating with the star. In 
this frame, the disk is being forced to push through and 
break all stellar held lines each rotation period. 

The result is that as the star spins, the held loops are 
forced to bunch up on one face of the disk and become 
rarihed on the opposed side of the disk (red diagram 
in Figure 2), similarly to how water is pressurized on 
the leading edge of a boat paddle. Unlike the azimuthal 
induction, this paddle-like braking occurs over the entire 
disk as opposed to solely the region orbiting beyond co¬ 
rotation. As such, this effect constitutes a signihcant 
source of angular momentum loss for the star. 

In addition to the intuitive braking torque upon stellar 
spin rate, the vertical component of magnetic induction 
leads to an additional torque affecting stellar orientation. 
In order to calculate the magnitude of these torques, we 
take a similar approach to that used within the azimuthal 
framework above. Specihcally, we suppose that the com¬ 
ponent of the magnetosphere parallel to the disk is built 
up over half the disk and reduced on the other half to a 
degree which is of the same order as the stellar magne¬ 
tosphere. Again, we suppose non-steady reconnection to 
provide a bound on magnetic field magnitude. Thus, we 
suppose the held induced by vertical held dragging (with 
subscript “p”) is given by 

= C{4') \Br Gr + , (26) 

where C(<^) = sin((^ — -0). Recalling that ip is the az¬ 
imuthal angle of the stellar spin pole direction projected 
onto the disk plane, this functional form for 0 ensures 
that held lines bunch up on the faces of the disk where 
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magnetic field lines are being pushed into the disk and 
the field is rarified on the other faces, where the field is 
being rotated away from the disk surface^. 

Finally, vertical field motion generates torques which 
take the form, 

7 7rT . , 

sm(/3) cos{'tp) 

TttT 

'^y = sm{tp) 

®in 

2ttT 

Tz = -w- 3 - cos(/3) sin(/3), 

(27) 

where we have again invoked Oout 3> Oin- 

We complete the specification of the torques arising 
from vertical field dragging by noting that if the star 
spins sufficiently slowly, the field lines will be able to 
diffuse vertically through the disk within one rotational 
period. Essentially, this is equivalent to saying that the 
field lines no longer reconnect and a true steady state is 
attainable (at least in terms of vertical field dragging). 
Suppose that the disk’s diffusivity is prescribed as 

r] = ash'^ riK, (28) 

where dg is the Shakura-Sunyaev parameter for disk vis¬ 
cosity (Shakura & Sunyaev 1973) and riK is the Keple- 
rian angular velocity. We may suppose that the stellar 
field lines could diffuse vertically through the disk over a 
timescale t ^ h'^/rj = l/(ds tt-k)- Put another way, if the 
star is spinning more slowly than about l/dg times the 
angular velocity of the gas in the disk, we would expect 
field lines to diffuse through with greatly reduced mag¬ 
netic induction. The likely value of dg within the inner 
disk ranges widely, from 10“^ to 10“^. Owing to this 
range, for definiteness we suppose that negligible field 
is induced if the star rotates less than about 10 times 
the Keplerian angular velocity of the disk. We include 
such an effect qualitatively by modulating the torques 
in equation (27) by a factor of exp[—(w/e Wq)”^] with 
e = 0.1 and uiq corresponding to an 8 day spin-period. 

2.4.3. Summary of Magnetic Torques 

The cumulative impact of the magnetic torques derived 
above is depicted visually in Figure 2 and may be sum¬ 
marized as follows. The inner regions of the disk (i.e. 
a' < a(.o) rotate faster than the stellar magnetosphere, 
dragging field lines azimuthally, ahead of the stellar ro¬ 
tation. The result is both an acceleration of the stellar 
spin rate and a torque that acts to realign the disk and 
stellar spin pole. Note that previous authors have con¬ 
cluded ambiguity over the sense of the tilting component 
to this torque (Lai et al. 2011) largely owing to uncer¬ 
tain parameterizations for the stellar wind. Our analysis 
here demonstrates that the pure impact of the magnetic 
torques from the inner disk is realignment. As can be 
seen from equations (25), the z-component (i.e., disk- 
aligned) component is always positive, which tends to 

^ The = sin(iji — ■0) form additionally ensures that V ■ Bp=0. 


spin up and realign. The horizontal components, how¬ 
ever are also positive and so tend to misalign, but simply 
adding the z- and horizontal components in quadrature 
results in a net realignment torque. 

The sense of torques is exactly reversed in the outer 
disk, i.e., the outer regions have a tendency to mis¬ 
align. Torques arising from vertical induction tend pre¬ 
dominantly to brake the stellar spin rate but they in¬ 
clude a small component acting to realign. In order 
to complete the specification of magnetic torques, we 
must choose a value for the inner disk radius a[^. In 
our models, we adopt a so-called “disk-locked” configu¬ 
ration (Konigl 1991; Mohanty & Shu 2008), whereby the 
outer and inner disk torques cancel at some prescribed 
stellar period (which we set to 8 days; see section (2.5) 
below). Accordingly, the inner disk truncation radius lies 
at ^ 0.8 Oco 1 / 3 = 0 ) in agreement with more sophisticated 
numerical simulations (Long et al. 2005). 


2.5. Stellar Spin Rate 

A complete account of the factors influencing the spin 
rates of young stars remains elusive (Matt & Pudritz 
2004; Herbst et al. 2007; Littlefair et al. 2010). The rea¬ 
son is that if one naively assumes a T-Tauri star to con¬ 
tract along its Hyashi track whilst conserving angular 
momentum, it would be expected to spin up sufficiently 
quickly such as to reach break-up velocity within the disk 
lifetime. However, observations clearly demonstrate that 
the vast majority of T-Tauri stars spin at rates far re¬ 
duced from their break-up angular velocity. Specifically, 
most young stellar objects rotate at about 3 — 10 day 
periods with only a rare few which are indeed close to 
break-up angular velocity, at periods of ^ 1 day. The 
slow spins of such stars is sometimes referred to as the 
‘stellar angular momentum problem’ because stars spin 
slower than they “should”. 

It has long been suspected that magnetic star-disk in¬ 
teractions play a key role in modulating the spin rates 
of T-Tauri stars (see Bouvier et al. (2014) for a recent 
review). Here, we do not attempt any sophisticated mod¬ 
eling aimed at explaining the stellar spin rates in detail. 
However, we note that the quasi-periodic twisting and 
opening of magnetic field lines discussed above provides 
a conduit through which angular momentum may be lost 
from the system, by way of stellar and/or disk winds. In 
addition, the torque arising from vertical field motion, 
derived above, provides an extremely efficient braking 
mechanism upon stellar spin, provided the stellar mag¬ 
netosphere is sufficiently strong. As such, mutual incli¬ 
nation between a star and its disk may give rise to hith¬ 
erto under appreciated magnetohydrodynamic influences 
upon stellar spin, especially when considering stars with 
lower masses (M < 1.2M 0 ) and strong magnetospheres 
(^1 kGauss). 

Within the framework of our model, we treat spin-rate 
evolution as follows. We allow the torques prescribed 
above to act freely on stellar spin rate. In addition, we 
introduce a relaxation factor, reflecting the reluctance 
of stars to spin up to break-up, whereby the star seeks 
the disk-locked spin rate Wr over a timescale Tr. The 
relaxation is prescribed as. 
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Fig. 2.— A schematic to illustrate the origin of each magnetic torque. The blue region represents the disk material interior to corotation, 
including the inner wall of the disk, which super-rotates with respect to the stellar spin, acting both to spin the star up and realign its axis 
with the disk’s. The green region is further out and so rotates more slowly, braking the stellar spin and acting to misalign the star. Red 
lines represent the requirement that the entire magnetosphere must be dragged vertically through the disk once per stellar rotation period 
if the star and disk are misaligned. A simple illustration of the physical mechanism behind each torque is shown on the right. The colored 
arrows in the top-left denote the net torque acting upon the stellar spin axis: green regions slow down and misalign the star, blue regions 
speed the star up and force realignment, whereas red regions act to brake stellar rotation whilst realigning the stellar spin-axis. Summed 
together, the resultant magnetic effect is usually to realign the disk and star with each other. 


doj 

dt 


UJ — UlQ 

5 

Tr 


(29) 


where we choose, for the stellar contraction timescale, 
which may be derived from equation (2): 


T, = [R/R)-^ 


24 ^ 

5-C GM^ 


(30) 


The above timescale changes as the star contracts, but 
remains of order ~ 1 Myr which is the value we adopt 
in this work. As already stated above, we choose the 
equilibrium angular velocity to be within the observed 
range, Wq = 27r/(8days). 


while assuming exclusively principal axis rotation (this is 
an excellent approximation for a T-Tauri star, spinning 
at a period of 3 — 10days). Specifically, jd is the angle 
between the central star’s spin axis and the binary orbit 
normal; Cl is the longitude of ascending node of the star 
in the binary frame where = 0 implies commensurate 
disk and stellar lines of nodes; and the third Euler angle 
(fi is the angle through which the star rotates as it spins 
{(fi only enters the equations as a rate of change: ip = uj). 

The equations for the evolution of /3, w and fl, adapted 
from Peale et al. (2014) are: 


2.6. Frame of Reference 

It is beneficial to carry out all calculations in the frame 
of a distant, binary companion to the central star. As 
such, we follow the approach of Peale et al. (2014) and 
define Euler angles within the binary frame related to 
the nutation, precession and rotation of the rigid body 
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duj 

dt 


0 

dt 


dCt 

dt 


W — Wo 


— Ny cos(ii) 


sin(;9) 


Nx sin(fi) 


■ cos(/3) Ng 


cos{i3){—Nx sin(il) + Ny cos(ii)) 


+ NgSm{P)) 

1 


sin(/3) 


Nx cos(fi) + Ny sin(ii) 


(31) 


where -/V^ are projected torques. Note that by fixing the 
disk’s longitude of ascending node at fl' = 0, we have im¬ 
plicitly placed ourselves into a frame co-precessing with 
the disk’s angular momentum vector, as discussed above. 
The effect of precession was included within the gravita¬ 
tional part of the calculation (equation (10)) and so we 
need not retain it here. 

The projected quantities are directly related to the 
torques calculated above, although the components of 
the torques in the disk frame, —r^/, must first be pro¬ 
jected onto the Cartesian axes in the binary frame. Such 
a projection constitutes a simple geometric rotation of 
co-ordinates because the disk-binary inclination is a con¬ 
stant of motion. The co-ordinate rotation transformation 
angle is fixed at some prescribed value, /3' to the z-axis 
and its sense is defined as anti-clockwise about the x- 
axis. As such, the components, N^ are given in terms of 
Ti by: 


Nx = -Tx'/{IM^RI), 

Ny = -(cos(/3')v - sin(,5'))/(/M*i?J), 

Nx = -{cos{I3')tx' + V sin(/3'))/(/M*i?J). (32) 

The above equations can be used to model the dynamics 
of the central star resulting from its magnetic interac¬ 
tions with its protoplanetary disk whilst the disk itself 
precesses within the binary frame. 

3. RESULTS 

We motivate the following analysis by considering the 
timescales over which the magnetic torques act, as de¬ 
rived above. Specifically, magnetic torques increase as 
oc meaning that the order-of-magnitude difference in 
field strengths between high and low-mass stars trans¬ 
lates to an enhancement of two orders of magnitude 
in the torques felt by lower-mass stars. Supposing the 
star to be set up in a retrograde, anti-aligned state 
(/? = tt), we calculate the characteristic realignment 
timescale Taiign using the equations derived above: 

^ _w 3GM^I /io 

J align - ^ ^ 2^4 (33 j 

where i?* follows the time-dependence described in equa¬ 
tion 3, giving rise to a time dependent magnetic torquing 
timescale which we illustrate in Figure 3. Taking 
1 kGauss as the held strength typical of low-mass stars, 


under nominal star-disk parameters their absolute re¬ 
alignment timescales are of the order of ^ 1 Myr through¬ 
out the majority of the disk lifetime. Conversely, the 
analogous timescale for the 0.1 kGauss helds typical of 
high-mass stars is closer to ^ 100Myr (Figure 3). These 
timescales are respectively shorter and longer than the 
typically-quoted 3 — 10 Myr lifetimes of protoplanetary 
disks (Haisch et ah 2001). Therefore, magnetic inter¬ 
actions only have the potential to wipe out primordial 
star-disk misalignments in low-mass systems. 



Fig. 3.— The approximate magnetic torquing timescale (Taiign) 
as a function of disk-star age for four regimes. The green lines ap¬ 
ply to high-mass stars with a surface field strength of ~ 0.1 kGauss. 
Red denotes low-mass stars with a surface field strength of ~ 
1 kGuass. In both cases, the upper line considers the timescale 
relevant to a star which is spinning with a 10 day period whereas 
the lower line applies to one with a 3 day period. The stellar spin 
rate is assumed to result from the locking to a circumstellar disk 
(Konigl 1991) and so the faster-spinning cases consider a disk which 
is truncated at smaller radii, increasing the magnetic influence. No¬ 
tice that the dominant effect upon magnetic torquing timescale is 
the magnetic field strength, with timescales proportional to its in¬ 
verse square. The timescales increase with time because the star 
contracts, leading to an effectively weaker field at the position of 
the inner disk. Only for the very earliest stages of protoplanetary 
disk evolution are high-mass stars’ magnetospheres strong enough 
to significantly alter their orientation whereas low-mass stars re¬ 
main dynamically influenced by magnetic fields throughout the en¬ 
tire typical disk lifetime of ~ 1 — 10 Myr. 


3.1. A companion inclined by 30 degrees 

In order to place magnetic realignment into the disk- 
torquing context, we integrate the mutual star-disk incli¬ 
nation over a 10 Myr timescale (i.e. an upper bound on 
typical disk lifetimes; Haisch et al. 2001), adopting the 
case of a 1 Mq massive companion inclined at 30 degrees 
to the disk plane. We present three cases in Figure 4: 
the purely-gravitational scenario (thick, pink line), the 
dynamics of a star with a dipole magnetic field strength 
of 0.1 kG, characteristic of higher-mass stars (thin grey 
line) and the situation for a dipole field strength of 1 kG, 
typical of low-mass stars (blue line). 

Both the magnetic field-free and 0.1 kG cases look al¬ 
most identical. In other words, the magnetic fields of 
higher-mass stars are dynamically unimportant within 
the disk lifetime, in accordance with the timescale anal¬ 
ysis quoted above. The gravitational evolution is then 
equivalent to that discussed in Spalding & Batygin 
(2014). Specifically, at the earliest times, the star pre- 
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cesses about the disk plane much faster than the disk pro¬ 
cesses about the binary plane (the high-frequency wig¬ 
gles at time t < 0.3 Myr) and small misalignments are 
maintained. However, as the star shrinks and the disk 
loses mass, gravitational star-disk coupling weakens until 
the two precession timescales are roughly commensurate. 
This situation causes the system to pass through a secu¬ 
lar resonance (Murray & Dermott 2000; Morbidelli 2002), 
facilitating a brief period of extremely efficient angular 
momentum transfer between the disk and star, resulting 
in large misalignments. 

Importantly, even with the initial companion-disk in¬ 
clination set to 30 degrees, a retrograde disk may emerge 
from the secular resonance. It is this non-linearity be¬ 
tween disk-binary inclination and resulting misalignment 
which is the essence of how disk-torquing can account for 
the entire range of observed misalignments (Spalding & 
Batygin 2014). 

The same secular resonance is encountered in the 1 kG 
case, but crucially, the resulting large misalignments are 
erased within typical disk lifetimes. An important as¬ 
pect of the dynamics is that the orientation of the star 
does not converge upon the disk plane, but rather, it con¬ 
verges to 30 degrees, i.e., the binary plane., which may be 
explained qualitatively as follows. Gravitational systems 
are conservative, allowing their dynamics to be described 
as following contours of a scalar function, or Hamiltonian 
(Morbidelli 2002). Magnetic torques introduce a dissipa¬ 
tive component that acts to turn the elliptical equilib¬ 
rium points of the Hamiltonian into attractors. In this 
case, the attractor is the binary-aligned state. 

A crucial point to emphasize is that in aligning with 
the binary plane, retrograde disks are prohibited. To see 
this, note that the companion’s gravitational perturba¬ 
tion upon the disk is dynamically equivalent whether the 
binary orbit is clockwise or counter-clockwise (as defined 
in some arbitrary frame). Such symmetry arises from ap¬ 
proximating the orbit as a massive wire and has the con¬ 
sequence that the equilibrium of the Hamiltonian never 
lies at a position where the star is tilted by more than 
90 deg to the disk. 

3.1.1. Sensitivity to Disk-Binary Inclination 

For completeness, we investigate the star-disk dynam¬ 
ics over a range of disk-binary inclinations in addition to 
the 30 deg case above. Specifically, Figure 5 illustrates 
the more extreme cases of 45, 60 and 75 degree incli¬ 
nations. The picture is very similar across 30, 45 and 
60 degrees, i.e., the magnetic fields of low-mass stars are 
capable of realigning them with the binary plane within 
the disk lifetime for a broad range of angles. In contrast, 
even the stronger fields were unable to wipe out retro¬ 
grade disks when an extreme initial binary-disk inclina¬ 
tion of 75 degrees is chosen. As expected, the orientation 
of the star is almost entirely unaltered by magnetic disk- 
star interactions for 0.1 kG fields in all cases. 

We also present the time evolution of stellar spin result¬ 
ing from the dynamics. A notable effect of larger binary 
inclinations is the considerable impact upon stellar spin 
rate, displayed in Figure (6). Though one should not 
take the rotation periods displayed in Figure (6) too lit¬ 
erally, owing to the uncertainties in prescribing rotation 
rates, an important aspect of the set-up emerges. Specif¬ 
ically, tilted stars will interact with their disks in such a 


way as to significantly brake stellar rotation, with such 
braking being more significant for larger star-disk incli¬ 
nations. The origin of such an effect is in the requirement 
of a tilted star to drag its entire magnetosphere vertically 
through the disk once every period^. For high binary in¬ 
clinations, the star is potentially spun down to nearly a 
stand-still as the system passes through the secular reso¬ 
nance, before being spun-up again in the direction of the 
disk. 

4. DISCUSSION 

4.1. Summary 

Prior to this work, a deficiency of the “disk-torquing” 
model (Batygin 2012; Batygin & Adams 2013; Spald¬ 
ing & Batygin 2014) for the acquisition of spin-orbit 
misalignments lay in its inability to reproduce the ob¬ 
served mass-dependence of misalignments (Winn et al. 
2010; Albrecht et al. 2012). Here, we resolved this dis¬ 
crepancy through the addition of a comprehensive set 
of magnetic star-disk torques on top of the previously- 
considered gravitational dynamics. Taking account of 
the recently-observed mass-dependence of T-Tauri dipole 
field strengths (Gregory et al. 2012), the observed trend 
arises naturally. 

We began by deriving analytical expressions describing 
the magnetic star-disk torques. Many previous authors 
have contributed to the development of a description of 
such torques (Ghosh & Lamb 1978; Livio & Pringle 1992; 
Armitage & Clarke 1996; Agapitou & Papaloizou 2000; 
Uzdensky et al. 2002; Matt & Pudritz 2004; Lai et al. 
2011). A common limitation has been that torques are 
considered to arise solely as a result of magnetic induc¬ 
tion associated with relative azimuthal motion between 
the stellar magnetosphere and the disk fluid. However, 
here also considered magnetic induction arising from 
oblique rotation. To this end, we have found that the 
associated torques provide an important correction to 
the existing pictures. Accordingly, this effect constitutes 
an efficient source of braking upon the stellar spin. 

The largest contribution to the torques affecting stellar 
orientation, as opposed to spin rate, arise from azimuthal 
induction. The influence of the disk material inside coro¬ 
tation is to spin the star up and align the stellar spin pole 
with the disk, with the opposite effect arising from the 
outer disk. This effect is amplified by stellar obliquity. 
In other words, when the star is tilted, a greater propor¬ 
tion of the disk is acting to accelerate the stellar spin and 
align it with the disk. Inclusion of this aspect increases 
the stability of the disk-locked equilibrium proposed else¬ 
where (Konigl 1991; Mohanty & Shu 2008). 

The final addition to our picture of magnetic torques 
was to include (to leading order) the effect of finite disk 
thickness. Whilst simple and widely adopted, a short¬ 
coming of the razor-thin disk model is that a star in¬ 
clined by 90 degrees to the disk is predicted to feel no 
magnetic torques arising from azimuthal field dragging 
alone (vertical induction still occurs). This is an artifact 
that comes about by supposing that there exists a negligi¬ 
ble solid angle within which field lines penetrate radially 
the disk. However, though disk aspect ratios are small 

^ Note that ignoring the torques arising from vertical field mo¬ 
tion would actually predict spin-up of a tilted star owing to the 
expanded corotation radius. 
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Fig. 4.— The time evolution of mutual star-disk inclination. We consider a binary companion to orbit the system at an inclination of 
30 degrees relative the disk (greater angles are displayed in Figure 5). The companion is prescribed to cause the disk to process with a 
1 Myr period. The purely gravitational case is shown as a thick, pink line. The thin, black line denotes evolution in the presence of the 
weak fields of high-mass stars (~0.1kGauss) and the blue line denotes the evolution corresponding to the strong fields of low-mass stars 
1 kGauss). As in previous work (Spalding & Batygin 2014), we find that a secular spin-orbit resonance is encountered as the disk loses 
mass and the star contracts. However, significant misalignments are inhibited for the stronger magnetic fields characteristic of low-mass 
stars whereas the fields of high-mass stars make no appreciable difference to the dynamics. Interestingly, the action of magnetic torques 
takes the star into alignment with the binary plane, not the disk plane, suggesting that small misalignments are indeed a natural outcome 
of low-mass star evolution, whereas high-mass stars can take on the full range of misalignments, in direct agreement with the observations. 


(~ 0.1), the inner edge of the disk is the closest point 
on the disk to the star and therefore can communicate a 
non-negligible torque to the central star (Equation 25). 

We motivated the potential for magnetic torques to 
sculpt the mass-dependence of obliquities by show¬ 
ing that the calculated magnetic torquing timescale is 
^IMyr for the IkG fields typical of low-mass stars 
whereas the analogous timescale applicable to high-mass 
stars was closer to 100 Myr. The fact that a typical pro¬ 
toplanetary disk lifetime sits right in the middle of these 
two timescales, at 3-10 Myr, indicates qualitatively differ¬ 
ent evolutionary scenarios during the disk-hosting stages 
of high and low-mass stars. In other words, star-disk 
magnetic torques constitute an influential factor in low- 
mass systems but may be neglected in high-mass sys¬ 
tems. 

The timescale analysis is compelling, and demonstrates 
that, were one to place a star and its disk in a misaligned 
configuration and leave them alone, a 1 kG dipole field 
would realign them before the disk dissipates. However, 
stars are likely to be aligned with their disks in the ab¬ 
sence of outside influences (Spalding et al. 2014). The ex¬ 
citation of misalignments must then come about by way 
of perturbations arising from a companion star (Baty¬ 
gin 2012). Appropriately, next we coupled the updated 
magnetic torques derived here to the purely gravitational 


model presented in Spalding & Batygin (2014). Specifi¬ 
cally, the star and the disk interacted through both gravi¬ 
tational and magnetic torques, with a massive companion 
causing the disk to process with a period of 1 Myr. 

As expected from the timescale analysis, the 0.1 kG 
fields characteristic of high-mass stars resulted in lit¬ 
tle deviation from a purely gravitational picture. How¬ 
ever, the stronger IkG fields caused the star to align, 
not with the disk, but with the binary plane (except for 
extreme disk-binary inclinations of >75 degrees). Thus, 
the most important result of our paper emerges: Dipole 
fields typical of lower-mass T Tauri stars realign the stel¬ 
lar spin axis with the plane of a perturbing companion on 
a timescale shorter than the disk lifetime. In this way, 
small misalignments occur naturally, but retrograde disks 
are prohibited in all but the more massive systems. 

Our results agree with the observational data presented 
in Figure 1. Retrograde circular orbits only exist among 
the higher-mass (M >1.2M0) population of stars. De¬ 
spite the absence of retrograde orbits around the lower- 
mass stars, significant obliquities still persist in general. 
It is important to note, also, that the observations only 
reveal the projected obliquity of exoplanetary systems. A 
true inclination of 60 degrees, say, will on average appear 
sky-projected as a lower obliquity within our current ob¬ 
servational dataset. On the other hand, distinguishing 
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Fig. 5.— Star-disk misalignments as functions of time for a variety of disk-binary inclinations. This set of figures show similar information 
to that shown in Figure 4, except we now illustrate the evolution for a range of angles: 30deg (top-left), 45deg (top-right), 60deg (bottom- 
left) and 75deg (bottom-right). The thick, pink line is the field-free case, the thin, black line presents the weak-field (O.lkGauss) case 
inherent to high-mass stars and the blue line denotes the strong-field (IkGuass) case of low-mass stars. The over-all pattern is largely 
similar up to 60 deg, in that the star is drawn towards a binary-aligned state over the magnetic realignment timescale. Not even the strong 
fields can undo the extreme resonant acquisition of misalignments occurring as a result of a 75 deg binary inclination. These larger binary 
inclinations are less likely, but raise the possibility that we may find a rare population of retrograde planetary orbits around low-mass stars 
in future datasets. 


retrograde from prograde orbits is a much simpler task. 
The lack of retrograde systems, together with the de¬ 
generacy inherent to measuring sky-projections, means 
that the current observations are fully consistent with 
our analysis here whereby small alignments persist while 
retrograde orbits are prohibited. 

A more subtle point is that, although the low-mass in¬ 
clinations approach binary-aligned equilibria, their time- 
evolution traces an oscillatory trajectory. Such oscil¬ 
lations in stellar orientation, especially for the 60 and 
75 degree cases (Figure 5), take the star retrograde with 
respect to the disk for a small portion of the trajectory. If 
the disk were to rapidly dissipate before the amplitude of 
oscillations was sufficiently damped by magnetic interac¬ 
tions, and at a point in the oscillations where the star was 
in a retrograde position, such a system could conceivably 
produce retrograde planetary systems around low-mass 
stars. This scenario is much less likely than the alter¬ 
native result of prograde systems, but does allow for the 
future detection of retrograde systems around low-mass 
stars, albeit at a significantly reduced frequency as com¬ 
pared to prograde systems. 

Our analysis was largely based upon order-of- 


magnitude estimates. In accordance with this level of 
precision, we had to make several assumptions about the 
star-disk configuration. Perhaps most importantly, we 
chose the inner truncation radius of the disk such as to 
mimic a “disk-locked” scenario whereby the inner and 
outer disk torques cancelled at a specified equilibrium, 
aligned spin rate. In accordance with observations (Bou- 
vier 2013), we decided upon 8 days for this equilibrium 
spin period in our simulations. As can be seen in Fig¬ 
ure (3), the chosen equilibrium period does indeed influ¬ 
ence the characteristic magnetic torquing timescale. Had 
we chosen 3 days, for example, disk-locking would require 
a smaller truncation radius, resulting in stronger disk- 
star torques (both magnetic and gravitational). How¬ 
ever, the relative variation in timescales between 3 and 
10 day configurations is far smaller than the difference 
between high and low-mass stars (owing to the de¬ 
pendence) . Furthermore, because 8 days is in the middle 
of these extremes, the end result with respect to star- 
disk inclination is qualitatively similar across all realistic 
star-disk configurations. 

One final omission that is worth a brief mention is 
the potential argument that, because our picture (disk- 
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Fig. 6.— The time evolution of absolute stellar angular velocity 
in the case of a 1 kGauss dipole field, plotted in units of the equi¬ 
librium angular velocity (here corresponding to an 8day period). 
We show the evolution for four values of disk-binary inclination, 
increasing from top to bottom: 30deg (cyan line), 45deg (green 
line), 60deg (purple line), 75deg (orange line). In each case, the 
star is prescribed to relax to a disk-locked equilibrium at an 8 day 
rotational period over a Kelvin-Helmholtz time. The relaxation is 
included in an ad hoc fashion and so the exact form of the rotation 
curves after ~ 2 Myr is not to be taken too literally. However, as is 
most apparent for the higher inclinations, magnetic braking con¬ 
stitutes a significant mechanism for the removal of stellar angular 
momentum. Indeed, for large binary inclinations, the star can be 
almost entirely stopped and re-spun within a relatively brief time 
interval. 

torquing) requires the presence of a companion, perhaps 
the mass-misalignment trend is merely a reflection of 
higher binarity in higher-mass stars. Whilst this particu¬ 
lar aspect has not been examined in depth, Crida & Baty¬ 
gin (2014) analyzed the expected spin-orbit distribution 
arising out of the purely-gravitational disk-torquing pic¬ 
ture presented in Batygin (2012). Whilst the expected 
distribution was consistent with that observed it was dif¬ 
ficult to simultaneously fit the highest and lowest obliq¬ 
uities and little attention was paid to any mass depen¬ 
dence. It is difficult to see how such a sharp transition in 
obliquities may arise purely from such a smooth occur¬ 
rence relation as mass vs binarity and so we feel the much 
sharper nature of the dipole field strength transition is a 
more likely cause. 

4.2. Viability of the Disk-Torquing framework 

Ultimately, our work here is a demonstration that the 
observed spin-orbit misalignments in exoplanetary sys¬ 
tems can come about predominantly by way of mech¬ 
anisms occurring during the early, disk-hosting phase. 
We do not expect the entirety of the observed misalign¬ 
ments to have originated through disk-torquing because 
dynamical interactions probably play a secondary role in 
the hot Jupiter delivery process (Ford & Rasio 2006). 
It is important, therefore, to contrast the predictions of 
each model in order to determine which mechanism, if 
any, is dominant. 

Dynamically-excited inclinations are likely to coincide 
with high eccentricities. Although controversy exists be¬ 
tween analytic and numerical analyzes over the sense of 
eccentricity evolution during disk-driven migration (e.g. 
Goldreich & Sari 2003; Bitsch & Kley 2010; Kley & Nel¬ 
son 2012), the general consensus is that eccentricity is 

® Indeed, general relativity provides a considerable boost to sta- 


limited to small values. As such, we predict that as 
datasets become more complete, the mass-misalignment 
trend will remain for circular systems, but not for eccen¬ 
tric systems (e > 0.1). Such a pattern is already begin¬ 
ning to emerge (Figure 1) but will be fully tested when 
the upcoming TESS mission commences. Furthermore, a 
disk-torquing origin for misalignments is fully consistent 
with the existence of multi-transiting systems (Huber et 
al. 2013). Indeed, we expect future observations of such 
systems to reveal the same mass-misalignment trend as 
that currently measured for single planet systems, with 
the break between high and low obliquities occurring at a 
similar (although not necessarily identical) stellar mass. 

Our work here has focussed largely on close in systems. 
What can be said about more distant planets? An im¬ 
plicit assumption we have made is that throughout the 
dynamics, the disk acts like a rigid, planar body. Nu¬ 
merical simulations have demonstrated that in reality, 
the disk is likely to develop a mild warp in response to 
the perturbations of a companion (Larwood et al. 1996). 
Therefore, the possibility exists for the outer planets to 
occupy a different plane from the inner planets despite 
having undergone disk-torquing. 

A more subtle difference between the outer and inner 
regions exists. One might imagine that close-in plan¬ 
ets, having migrated through their natal disks, will be 
much more shielded from dynamical instabilities arising 
from perturbations later in the systems lifetime. This 
shielding arises both as a result of sitting in a deeper 
gravitational potential well, but also through the stabi¬ 
lizing influence of general relativity®. Accordingly, we 
might expect there to exist a population of close-in sys¬ 
tems which have obtained their orbits by way of copla- 
nar disk migration, in addition to a separate population 
of planets around the same star that, for whatever rea¬ 
son, remained at larger radii within their disk, becoming 
subject to various dynamical instabilities later in their 
evolution. 

In conclusion, the magnetically-facilitated realignment 
presented here provides a natural pathway for the gen¬ 
eration, not only of high obliquities, but also their ob¬ 
served dependence upon mass. Such a “disk-torquing” 
model appeals to no assumptions beyond those within 
the bounds of what has been observed for T-Tauri sys¬ 
tems, namely magnetic field strength, stellar spin rate, 
disk ionization state and disk mass. Out of such nom¬ 
inal parameters, the observed mass-misalignment trend 
arises naturally. Additionally, the conclusions presented 
here are fully consistent with the decades-old picture of 
smooth migration through a protoplanetary disk (Gol¬ 
dreich & Tremaine 1980), leading to short-period orbits 
coplanar with the disk. 
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